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Abstract 

We study the properties of non-rotating and rotating neutron stars for a new set of equations 
of state (EOSs) with different high density behaviour obtained using the extended field theoretical 
model. The high density behaviour for these EOSs are varied by varying the ui— meson self-coupling 
and hyperon-meson couplings in such a way that the quality of fit to the bulk nuclear observables, 
nuclear matter incompressibility coefficient and hyperon-nucleon potential depths remain practi- 
cally unaffected. We find that the largest value for maximum mass for the non-rotating neutron 
star is 2.1M & . The radius for the neutron star with canonical mass is 12.8 — 14.1 km provided only 
those EOSs are considered for which maximum mass is larger than 1.6M@ as it is the lower bound 
on the maximum mass measured so far. Our results for the very recently discovered fastest rotating 
neutron star indicate that this star is supra massive with mass 1.7 — 2.7 Mq and circumferential 
equatorial radius 12 — 19 km. 

PACS numbers: 26.60. +c,91.60.Fe,97.10.Kc,97.10.Nf,97.10.Pg 
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I. INTRODUCTION 



The knowledge of neutron star properties is necessary to probe the high density behaviour 
of the equation of state (EOS) for the baryonic matter in (3— equilibrium. The EOS for the 
densities higher than p — 0.16 fm~ 3 can be well constrained if radii for the neutron stars 
over a wide range of their masses are appropriately known. Even the accurate information on 
the maximum neutron star mass M max and radius -R1.4 for the neutron star with canonical 
mass (lAMp) would narrow down the choices for the plausible EOSs to just a few. Till 

nnn 

date, the neutron stars with masses only around 1.4M are accurately measured |2J, 13J. 
Recent measurement of mass of the pulsar PSR J0751+1807 imposes lower bounds on the 
maximum mass of the neutron star to be 1.6M and 1.9M with 95% and 68% confidence 
limits, respectively [4|. The increase in the lower bounds of the neutron star maximum mass 
could eliminate the family of EOSs in which exotica appear and substantial softening begins 
around 2 to 4 p leading to appreciable reduction of the maximum mass. The available 
data on the neutron star radius have large uncertainties . The main source of 

the uncertainties in the measurements of the neutron star radii are the unknown chemical 
composition of the atmosphere, inaccuracies in the star's distance and high magnetic fielc 
(~ 10 12 G). Recent discovery of the binary neutron star system PSR J0737-3039A,B [2] 
with masses of the individual star being 1.338M and 1.249M have raised the hope for the 
possibility of measuring the moment of inertia due to the spin-orbit coupling effects [3j. It 
is expected that a reasonably accurate value for neutron star radius can be deduced from 
the moment of inertia measurements. Very recent discovery of the fastest rotating neutron 
star with rotational frequency of 1122 Hz observed in the X-ray transient XTE J1739-285 
[lo| has placed an additional constrained on the EOS at very high density 11]. 

Theoretically, the mass-radius relationship and compositions of the neutron stars are 
studied using various models which can be broadly grouped into (i) non-relativistic potential 
models 12] , (ii) non-relativistic mean-field models 



relativistic mean-field models (FTRMF) 



model 
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16] , (iii) field theoretical based 



19| and (iv) Dirac-Brueckner-Hartree-Fock 
231 ] . Each of these models can yield EOSs with different high density 



behaviour which is not yet well constrained. As a result, neutron star properties vary over 
a wide range even for the same model. In this work we shall mainly focus on the variations 
in the properties of the neutron stars obtained within the FTRMF models. The FTRMF 
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models predict the values of M, 



neutron stars 



3, 
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= 1.2 — 3.0Mq and R\a = 10 — 16 km for the non-rotating 
26j. The lower values of M max and R\a correspond to the neutron stars 



composed of nucleons and hyperons in {3— equilibrium, where as, the higher values of M max 
and i?i.4 correspond to the neutron stars with no hyperons. We would like to emphasize that 
not all the different parameterizations of the FTRMF model, employed to study the neutron 
star properties, are able to reproduce satisfactorily the basic properties of finite nuclei and 
nuclear matter at the saturation density. For instance, the value of the nuclear matter 
incompressibility coefficient which largely controls the low density behaviour of a EOS very 
in between 200 — 360 MeV for different FTRMF models. Though, the value of nuclear matter 
incompressibility coefficient is very well constrained to 230 ± 10 MeV by the experimental 
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The variations in 



data on the isoscalar giant monopole resonances in heavy nuclei 
the neutron star properties resulting from the differences in the high density behaviour of 
the different EOSs can be appropriately studied only if the low density behaviour for each of 
these EOSs are constrained using the experimental data on the bulk properties of the finite 
nuclei and nuclear matter at the saturation density. 



The extended FTRMF model 



29, 3q], 



31| includes mixed and self-coupling terms for the 



cr, uj and p mesons. The u;-meson self-coupling term enables one to vary the high density 
behaviour of the EOS without affecting nuclear matter properties at the saturation density 



24j. The mixed interaction terms involving p-mesons allow ones to significantly vary the 
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34j which plays crucial role 



density dependence of the symmetry energy coefficient 3^, 
in determining cooling mechanism of a neutron star [35]. Yet, such a versatile version of 
the FTRMF model is not fully explored to study the variations in the properties of the 
neutron stars resulting mainly from the uncertainties in the high density behaviour of EOS. 
In the present work we use the extended FTRMF models to obtain a new set of EOSs with 
different high density behaviour for the (3— equilibrated matter composed of nucleons and 
hyperons. Each of these different EOSs correspond to different choices for the uo— meson 
self-coupling and hyperon-meson couplings which mainly affects the high density behaviour 
of a EOS. The remaining parameters of the model are calibrated using a set of experimental 
data on the total binding energy and charge rms radii for a few closed shell nuclei. In our 
calibrational procedure we also use the value of neutron-skin thickness for the 208 Pb nucleus 
as one of the data. Since, the neutron-skin thickness is only poorly known, we obtain different 
parameter sets for different neutron-skin thickness ranging from 0.16 — 0.28 fm. We further 
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restrict the parameters to yield a reasonable value for the nuclear matter incompressibility 
coefficient at the saturation density. We use our EOSs to study the mass-radius relationship 
and chemical compositions for non-rotating neutron stars. For the case of rotating neutron 
stars, we present our results for the Keplerian sequences and also investigate the variations 
of mass and circumferential equatorial radius for the very recently discovered fastest rotating 
neutron star. 

In Sec. |TT]we outline very briefly the Lagrangian density and corresponding energy density 
for the extended FTRMF model. In Sec. II I II we present our various parameterizations for 
different combinations of a cu-meson self-coupling, hyperon-meson couplings and neutron- 
skin thickness for the 208 Pb nucleus. In Sec. II VI we present our results for the nuclear matter 
properties at the saturation density. In this section we also discuss about the quality of the 
fits to the finite nuclei for these parameterizations. In Sec. |V] we present our results for 
the properties of non rotating neutron stars. We also generate some rotating neutron star 
sequences for which the results are presented in Sec. IVII Finally our main conclusions are 
presented in Sec. IVHI 



II. EXTENDED FIELD THEORETICAL MODEL 



The effective Lagrangian density for the FTRMF model generally describes the interac- 
tions of the baryons via the exchange of a, u and p mesons. The a and the u mesons are 
responsible for nuclear binding while p meson is required to obtain the correct value for 
the empirical symmetry energy. The cubic and quartic terms for the self-interaction of the 
cr-meson are often considered which significantly improves the value of the nuclear matter 
incompressibility. Nevertheless, the value of the nuclear matter incompressibility coefficient 
for these models are usually larger in comparison to their values extracted from the exper- 
imental data on the isoscalar giant monopole resonances. Moreover, the symmetry energy 
coefficient and its density dependence is also somewhat higher relative to the corresponding 
empirical estimates. One can easily overcome these issues in the extended FTRMF model 
which includes self and mixed interaction terms for a, uj and p mesons upto the quartic 
order. In particular, mixed interaction terms involving rho-meson field enables one to vary 
the density dependence of the symmetry energy coefficient and the neutron skin thickness in 
heavy nuclei over a wide range without affecting the other properties of finite nuclei 33|, [34]. 
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The contribution from the self interaction of cu-mesons plays important role in varying the 
high den sity behaviour of the EOS and also prevents instabilities in the calculation of the 



EOS 
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361 ] . On the other hand expectation value of the p-meson field is order of magnitude 



smaller than that for the cu-meson field 



3l| . Thus, inclusion of the p-meson self interaction 



can affect the properties of the finite nuclei and neutron stars only very marginally [241 ] . 
The Lagrangian density for the extended FTRMF model can be written as, 

£ = C-BM + C a + Cu, + C p + Caup + C-em + £ep + LyY- (1) 

Where the baryonic and mesonic Lagrangian Cbm can De written, 

£bm = ^2^ B [il^d^ - {M B - g aB a) - {g uB l^n + ^9pBl^r B .p^ B . (2) 

B 

Here, the sum is taken over the complete baryon octet which consists of nucleons, A, E and 
H hyperons. For the calculation of finite nuclei properties only neutron and proton has been 
considered. r B are the isospin matrices. The Lagrangian describing self interactions for a, 
u>, and p mesons can be written as, 

C = \{d^a - m> 2 ) - |s> 3 - ^ N a\ (3) 
C» = + \ml^ + j£gt N {u^) 2 , (4) 



£p = ~\p^ lv + \m 2 pP ^ + ^gUip^f- (5) 



The uj^, p^ v are field tensors corresponding to the u and p mesons, and can be defined as 
^fiu _ qii^v _ Q v lo ^ anc [ pV-v — QVp v — d u p^. The mixed interactions of a, u, and p mesons 



£-au P can be written as, 



(6) 



+ ^3 91n9 2 pN^^P^ 
The C em is Lagrangian for electromagnetic interactions and can be expressed as, 

Cem = —FpuF^ ~ e^BjJ^p^A^B, (7) 

B 

where, F^ u = d^A 1 ' — d u A fl . The hyperon-hyperon interaction has been included by in- 
troducing two additional mesonic fields (ex* and <p) and the corresponding Lagrangian C Y y 



(Y = A, £, and H) can be written as, 



Y 



+ \ [d v a*d v a* - m>* 2 ) - l -S pv S^ + \m\^. 



(8) 



The charge neutral neutron star matter also includes leptons such as e~ and \x~ in addition 
to neutrons, proton, and hyperons at the high densities. The leptonic contributions to the 
total Lagrangian density can be written as, 



£=e,fi 



(9) 



The equation of motion for baryons, mesons and photons can be derived from the La- 
grangian density defined in Eq.fll]). The equation of motion for baryons can be given as, 

1 + T 3B 

<V - 94>B<P r J ~ 

(10) 



2 M 



{M B + g aB cr + g a *BV*) 



= e B ^ B . 



The Euler-Lagrange equations for the ground state expectation values of the mesons fields 



arc: 



(-A + ml ) a = 9<rBpsB ~ ^gl N ^ 2 - q^nC 3 + ®l9o 



n91n u2 



(11) 



+ ax&NglN™ 2 + a 29<rNg 2 pB p 2 + gl N g 2 pN (T p 2 , 



B 



(-A + ml) u = ^2 9ubPb ~ g^jv^ 3 - -<>j//<7 v«C.\ (7^ - <u'9a\ <C \ <r'^ 



-' „2 2 2 



(12) 



- a 3 'gl N g 2 pN^p 2 , 



(-A + m 2 p ) p = ^2gp B r 3 BpB - |^p7vP 3 - ^a 2 g aN g 2 pN ap - a 2 'gl N g 2 pN a 2 p 

B 

-r„2 „2 , ,2 



6' 



"3 9un9 p n^ P> 

(-A + m^») tr* = ^g^BpsB, 

B 

(-A + mJ)0 = J] 

B 



(13) 

(14) 
(15) 



- AA = ep p . 



(16) 
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where the baryon density pb, scalar density p S B and charge density p p are, respectively, 

6tt 2 ' 



PB = <*B7°*B> = ti > ( 17 ) 



^ = = ^/„ ^TPTW' (18) 

P P = • (19) 



\ 2 

Where, 7 is the spin degeneracy. The M B = Mb — g a B& — 9a* B a * is the effective mass of 
the baryon species B, ks is its Fermi momentum and t S b denotes the isospin projections of 
baryon B. 

The energy density of the uniform matter in the extended FTRMF models is given by 

1 f kj 1 1 

£ = ^ — k 2 ^jk 2 + Mfdk + 9uB^Pb + Yl 9 P BT 3B p + ^mla 2 
j=B,e 71 Jo B B 

1 K 3 3,^4 4 C4 4 £44 1 2 2 1 2 2 
+ + 24^^ - - -g pN p - -m^ - -m pP 

1 1 

- mg*N9* N (ruj 2 - -W'gl N gl N o- 2 u 2 - a^g aN g 2 pN ap 2 - -a^' gj N g 2 pN a 2 p 2 

~ \®z 91n9 2 p nU 2 P 2 + ^m 2 a ,a* 2 + Yd<t>B<f>PB - -^n\ft \ 

B 

The pressure of the uniform matter is given by 



l /C 4 rf/C loo K O O A 



A 3 - 2 Jo ^Mf 2 - n-' 2 1 



22 ^ 3 3 /v 4 4 



+ ^9t N ^ + ^pNP 4 + \ m ^ 2 + \ m lp 2 + ai9cNgl N <7u 2 ^ 



1 1 

' 2 2 _2, .2 , — „ 2 _ 2 , / 2 2 2 2 



+ 2 ftl 9aN9uN° U + a2g„Ng pN Op + ~«2 9 a N9pN^ P 
1 / 2 2 22 1 2 *2 . 1 2i2 

+ 9u>n9 p n^ P - -m a ,a + -m . 



III. PARAMETERIZATIONS OF THE EXTENDED FTRMF MODEL 

In this section we consider various parameterizations of the extended FTRMF model. The 
different parameter sets are obtained using different values for the w-meson self-coupling (, 
hyperon-meson couplings (i = a, u, p, a* and mesons) and neutron-skin thickness Ar 
for the 208 Pb nucleus. The parameter ( mainly affects the high density behaviour of the EOS 
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and can not be well constrained by the properties of the finite nuclei. The different sets of giy 
can be obtained to yield different EOSs for the dense matter without affecting the resulting 
potential depth for hyperons in the nuclear matter at the saturation density. The value of 
Ar for a single heavy nucleus like 208 Pb which can constrain the linear density dependence 
of the symmetry energy is only poorly known. The different choices for the (, giy and Ar 
are so made that they span entire range of values as often used in the literature. We must 
point out that the contributions from the p-meson self-coupling are ignored, because, their 
effects are found to be only marginal even for the pure neutron matter at very high densities 



24]. 



Towards our parameterizational procedure we first set hyperon-meson couplings g^y = 
in Eqs. (12151) . Then the remaining coupling parameters appearing in Eqs. (T5] - [6]) are 
determined by fitting the FTRMF results to the experimental data for the total binding 
energies and charge rms radii for many closed shell normal and exotic nuclei. We consider 
total binding energies for 16 ' 24 , 40 ' 48 Ca, 56 ' 78 Ni, 88 Sr, 90 Zr, 1 oo,n6,i32g n and 2o 8pb nudei; 
charge rms radii for 16 , 40,48 Ca, 56 Ni, 88 Sr, 90 Zr, 116 Sn and 208 Pb nuclei. In addition, we also 
fit the value of neutron-skin thickness for 208 Pb nucleus. Recently extracted value of neutron- 
skin thickness for 208 Pb nucleus from the isospin diffusion data lie within 0.16 — 0.24 fm 
indicating large uncertainties {37}]. We generate twenty one different parameter sets using 
different combinations of ( and Ar. The value of ( is taken to be 0.0, 0.03 and 0.06 and for 
the Ar we use 0.16, 0.18, 0.28 fm. The best fit parameters are obtained by minimizing 
the x 2 function given as, 



x 2 



1 Nd ( Q exp — Q th \ 2 



N d -N n ^\ S 



y p i=i 

where, Nd is the number of experimental data points and N p the number of parameters to 
be fitted. The Si stands for theoretical error and O exv and Of 1 are the experimental and 
the corresponding theoretical values, respectively, for a given observable. Since, the Of 1 in 
Eq. (T22T) is calculated using the FTRMF model, the value of x 2 depends on the values of the 
parameters appearing in Eq. (J2]-[n]). The theoretical error Si in Eq. (122|) are taken to be 1.0 
MeV for the total binding energies, 0.02 fm for the charge rms radii and 0.005 fm for the 
neutron-skin thickness. The best fit parameters for a given set of values of O exv and Si are 



searched using the simulated annealing method 38|, |39[]. In our earlier work 39] we have 



obtained the parameter sets for ( = 0.0,0.03 and 0.06 with Ar = 0.18 fm. Here too we 



follow the same strategy to obtain the parameter set for a given combination of Ar and £. 
In Tables [J [IT] and II III we list the values of parameters for all the sets presently generated. 

We now determine the values of the hyperon-meson coupling parameters g^y. These 
couplings can be expressed in terms of the nucleon-meson couplings using SU(6) model as, 

1 11 

-^9<tN — — ^crS — 9aE, 

1 11 

-^9wN — -j9uk — ^fl'wE — 9u>~, 

g P N = g P T, = %g P s, g P A = 0, 

2(7cr*A — 2(7 CT *s — ga*5 — g 9uiN, 9a*N — 0, 

2^a = = 9<f>3 = ~~^~9uNi 94>n = 0. (23) 

The neutron star properties are quite sensitive to the values of g„y and g w y. Where as 
neutron star properties do not get significantly affected even if the values of g a *y is varied 
over a reasonable range for a fixed value of g<py [40]. For g pY , g a * Y and g^y we use the 
values as given by Eq. fT23|) . The values of g a y and g^y are determined using the expressions 
for the hyperon-nucleon potential. The potential depth for a given hyperon species in the 
nuclear matter at the saturation density {p sa t) is given as, 

U^ipsat) = -9aY(?{psat) + 9cuY^(psat)- (24) 

The values of Uy^ chosen are as follows 41], 



U ( A N) = -28 MeV, U^ N) = +30 MeV and ui N) = -18 MeV. (25) 

Normally, g G y is determined for a given value of UyiPsat) with g^y taken from SU(6) 
model. For the sake of convenience we define, 



9mY 
9mN 



mY 



for A and X hyperons 

(26) 

2 ( 9mY_ \ f or ^ hyperons, 

^ \ 9mN J 

where, m stands for a and u mesons. In the present work we vary X w y from 0.5 — 0.8 
42J . In Fig. [1] we display the variations of X a y as a function of X^y obtained using the 
parameter set corresponding to ( = 0.03 and Ar = 0.22 fm. The values of X a y are for 
all other combinations of ( and Ar are very much the same as depicted in Fig. [TJ This is 



due to the fact that the properties of symmetric nuclear matter, like, binding energy per 
nucleon B/A, nuclear matter incompressibility coefficient K, effective nucleon mass at 
the saturation density are very much similar for all the parameterizations considered in the 
present work. 



IV. NUCLEAR MATTER AND FINITE NUCLEI 



The various properties associated with the nuclear matter are obtained using parameter 
sets of Tables HI HT1 and HTT1 The values of B/A, K, and p sat for all these parameter sets 
lie in a narrow range. We find that B/A = 16.11 ± 0.04 MeV, K = 230.24 ± 9.80 MeV, 
M* N /M N = 0.605 ± 0.004 and p sat = 0.148 ± 0.003 fm" 3 . The values of the symmetry energy 
coefficient J and its linear density dependence, 



T Q dJ 

L=3 % 



(27) 



Psat 

are strongly correlated with the Ar for the 208 Pb nucleus used in the fit. In Fig. [2]we display 
the variations of J and L calculated at saturation density as a function of Ar. The values 
of L lie in the range of 80 ± 20 MeV for Ar varying in between 0.16 to 0.28 fm which is in 



reasonable agreement with the recent predictions based on the isospin diffusion data [37]] . 

The relative errors in the total binding energy and charge rms radius for the nuclei 
included in the fits are more or less the same as we have obtained in our earlier work 39]. 



So, we do not wish to present here the detailed results. It might be sufficient for the present 
purpose to display the results for the rms errors for the total binding energies and charge 
rms radii obtained for our newly generated parameter sets. In Fig. [3] we plot the rms errors 
for the total binding energies and charge radii as a function of Ar. It is quite clear from 
this figure that rms error show hardly any variations implying that all the parameter sets 
generated in the present work fit the finite nuclear properties equally well. In fact, if we do 
not consider the parameterizations with ( = 0.0 and Ar = 0.26 or 0.28 fm, the rms errors 
on the total binding energy are 1.5 - 1.8 MeV which is comparable with one obtained using 



NL3 parameterizations as most commonly used 43j. The rms error of charge radii for the 



nuclei considered in the fit lie within the 0.025 — 0.040 fm. 
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V. NON-ROTATING NEUTRON STARS 



In this section we present our results for the properties of the non-rotating neutron stars 
for a set of EOSs obtained using different parameterizations for the extended FTRMF model. 
Each of these parameterizations corresponds to different combinations of neutron-skin thick- 
ness Ar in 208 Pb nucleus, the u -meson self-coupling ( and hyperon-meson couplings X^y as 
described in Sec. IIHI The values of Ar, ( and X u y are so varied that they span the entire 
range of values as often encountered in the literature. The variations in ( and X^y affect 
the high density behaviour of the EOS, whereas, the density dependence of the symmetry 
energy coefficient is strongly correlated with Ar. It is therefore natural to expect that the 
variations in Ar, ( and X w y can affect significantly the neutron star properties. The param- 
eters of FTRMF model are so calibrated that the quality of fit to finite nuclei, the properties 
of nuclear matter at saturation density and hyperon-nucleon potentials are almost the same 
for each of the parameterizations. Thus, these parameterizations provide the right starting 
point to study the actual variations in the properties of neutron star resulting from the 
uncertainties in the EOS of dense matter. 

The properties of non-rotating neutron stars are obtained by integrating the Tolman- 



Oppenheimer-Volkoff (TOV) equations [44(. To solve the TOV equations we use the EOS 
for the matter consisting of nucleons, hyperons and leptons. The composition of matter at 
fixed total baryon density, 



P = ^Pb, (28) 

B 



are so determined that charge neutrality condition, 

2?BPfl + 53?/P< = 0, (29) 

B t 

and chemical equilibrium conditions, 

HB = H>n- iB^e (30) 

Afy = At e (31) 

are satisfied. In Eqs. (1291131 \ q and fi are the charge and chemical potential for various 
baryons and leptons considered in our calculations. For densities higher than 0.5po, the 
baryonic part of EOS is evaluated within the FTRMF model. Whereas, the contributions of 
the electrons and muons to the EOS are evaluated within the Fermi gas approximation. At 

11 



densities lower than 0.5po down to 0.4 x 10 10 po we use the EOS of Baym-Pethick-Sutherland 
(BPS) |45|. 

In Fig. [4] we plot the EOS for the pure neutron matter and symmetric nuclear matter as a 
function of number density for the selected combinations of ( and Ar. We see that the EOS 
for ( = 0.0 is the most stiffest, and as ( increases the EOS becomes softer. The softening 
of EOS with ( is more pronounced at higher densities. In Fig. [5] we plot our results for 
the neutron and electron chemical potentials as a function of baryon density obtained for 
the EOSs corresponding to the moderate values of Ar and X w y The chemical potentials 
for other particles can be evaluated using Eqs. ( 1301) and ( I3T1) . The change in slope for 
neutron chemical potential vs. baryon density is associated with appearance of hyperons. 
The decrease in /i e for p ^ 2p is accompanied by the appearance of the S~ hyperons. The 
maximum values of fi e is less than half of the bare mass for kaons which indicate that the 
presence of hyperons inhibits the kaon condensation. 

Let us now consider various neutron star properties resulting from the EOSs for the two 
different parameter sets referred hereafter as LY and UY. These parameter sets are obtained 
using different combinations of Ar, ( and X^y- The parameters of LY set are obtained with 
Ar = 0.16 fm, ( = 0.06 and X^y — 0.5. Whereas, UY parameterization is obtained with 
Ar = 0.28 fm, £ = 0.0 and X^y = 0.8. Among all the parameterizations as obtained in 
Sec. IHIt LY and UY yield the softest and the stiffest EOSs, respectively. Thus, maximum 
variations in the neutron star properties can be studied using the EOSs obtained for LY 
and UY parameter sets. For the comparison, we also present our results for the L0 and U0 
parameter sets similar to LY and UY parameterizations, but, with no hyperons. In Fig. [6] we 
present our results for mass-radius relationship for LY, UY, L0 and U0 parameterizations. 
The region bounded by R ^ 3GM/c 2 is excluded by the causality limit [4g]. The line 
labeled by AI /I = 0.014 is radius limit estimated by Vela pulsar glitches j^]. The rotation 
constraint as indicated in Fig. [6] is obtained using [48[ , 

/ M \ 1/2 / 10km \ 3/2 TT , . 

^ 1833 HmJ hrj (32) 

with rj = 0.57 and Vf. = 1122 Hz which corresponds to the frequency for the fastest rotating 
neutron star present in the recently observed X-ray transient XTE J1739-285 [h]]. The 
renormalization factor rj account for the effects due to deformation and gravity. We also 
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calculate the variations in the radiation radius, 



R 



R 



(33) 



•oc 



2GM 
Re 2 



for the neutron star with the canonical mass 1.4M . It can be verified by using the results 
for the LY and UY cases presented in Fig. [6] that Roo lies in the range of 14.2 — 16.8 km. 
Similarly, without the inclusion of hyperons, the values of R^ vary in the range of 15.3 — 
16.8 km. In Tables IV and V we collect few important bulk properties for the non-rotating 
neutron stars with maximum and canonical masses. We see that the values of M max with the 
inclusion of hyperons varies between 1.4 — 2.1M . Once the contributions from the hyperons 
are ignored M max varies between 1.7 — 2.4M . The values of R\a varies from 11.3 — 14.1 km 
and 12.5 — 14.1 km depending on whether the hyperonic contributions are included or not. 
Thus, combining our results for the neutron stars with and without hyperons we find that 
the values of M max and R\a obtained within the FTRMF model can vary about 1M and 3 
km, respectively. These variations are almost half of the ones obtained earlier using FTRMF 
model in which bulk nuclear observables and nuclear matter incompressibility were not fitted 
appropriately. The values of redshift given in the Tables IV and V are obtained for the ratio 
M/R as, 



Our results for the values of redshift for the neutron star with canonical mass are 0.22 ±0.03. 
It is also interesting to note that Z ^ 0.35 only for the stars with masses 1.7M or larger. 

In Fig. [7] we have plotted the threshold densities for various hyperon species. In the same 
figure we also show the values of central densities for the neutron stars with canonical mass 
and maximum mass. The threshold density is lowest for the A hyperons. It is interesting to 
note that for the UY case the threshold density for the A hyperons is almost equal to the 
central density for the neutron star with the canonical mass. This implies that the properties 
of the neutron star with the canonical mass do not get affected by the hyperons for the UY 
parameterization. This is the reason that our results for the mass and radius for U0 and UY 
parameterizations are very much similar for the neutron stars with masses upto 1.6M as 
can be seen from Fig. [61 The S + and S° hyperons do not appear in density range relevant 
for the present study. However, for the TM1 parameterization of the FTRMF model one 



1 



(34) 



Z 




1 
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finds that all kinds of hyperons appear well below 7po |49|, |5Cj • This seems to be due to large 
value of nuclear matter incompressibility coefficient (K = 281 MeV) associated with the 
TM1 parameter set. In other words, not only the variations in the properties of the neutron 
stars reduces but also the chemical compositions for these stars can become different if the 
parameters of the FTRMF models are calibrated appropriately. In Fig. [S]we plot the particle 
fractions as a function of radial coordinate. These fractions are calculated for the neutron 
stars with M max = 1.4M Q and 2.1M corresponding to the LY (upper panel) and UY (lower 
panel) parameterizations, respectively. The neutron fractions in Fig. [8] are plotted after 
dividing them by a factor of three. We see that the compositions of the neutron stars shown 
in the upper and lower panels are not the same. For the case of LY parameterizations, H~ 
and £~ hyperons appear more or less simultaneously. For the UY case, S° hyperons appear 
instead of S~ hyperons . It is noteworthy that for the case with UY parameterization the 
hyperons are the dominant particles at the interior (r < 4 km) of the neutron star leading 
to complete deleptonization. We see from Fig. [8] that the proton fractions for both the cases 
are greater than the critical value (~ 15%) for the Direct Urea process to occur [351 ] . 

We now consider our results for neutron star properties at the canonical and maximum 
masses for the set of EOSs obtained using all the different parameterizations as given in Sec. 
IIHI These different parameterizations correspond to the different combinations of the Ar, 
( and X u y The values of Ar, ( and X^y vary in the range 0.16 — 0.28 fm, 0.0 — 0.06 and 
0.5 — 0.8, respectively. The knowledge of M max and (-R1.4) or moment of inertia (X L4 ) for 
neutron star with canonical mass are very important in order to understand the behaviour 
of the EOS over the wide range of density well above po- The discovery of the pulsars PSR 
J0737-3039A,B and PSR J075 1+1807 have raised hope for availability of more accurate 
information about these quantities in near future. The M max probes densest segment of the 
EOS. Whereas, R1.4 or X1.4 probes relatively lower density region of EOS. It is not possible 
to say a priori whether or not M max is correlated to the properties of neutron star with 
1.4M m . Earlier studies using FTRMF models indicate some correlations between M max and 
Ria [19]. Another study carried out for 25__EOSs taken from different models show hardly 



any correlations between M max and X x 34 



51] . In Fig. [9] we plot the variations of radius and 



the redshift for the neutron star with the canonical mass as a function of M max . We see 
that M max varies between 1.4 — 2.1M and -R1.4 varies between 11.3 — 14.1 km. The vertical 
line at M max = 1.6M corresponds to the mass of the PSR J0751+1807 measured with 95% 
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confidence limit. If only those EOSs are considered for which M max ^ 1.6M then the value 
of Ria would lie in the range of 12.8 — 14.1 km. This result is in reasonable agreement with 
R 1A = 14.8jl 1 ;| km as deduced very recently by adequately fitting the high quality X-ray 
spectrum from the neutron star X7 in the globular cluster 47 Tucanae [52]. We also note 
strong correlations of M max with R 1A and Z 14 . For a given value of M max , the spread in 
the values of Ria is 0.7 ± 0.1 km. Only for the M max ~ 1.4M Q we find that spread in the 
values of Ria is ~ 0.3 km. To understand it better, we list in Table VI the values of M max 
and Ria obtained for the sets of EOSs corresponding to the selected combinations of Ar, 
( and X^y For additional information we also give in Table VI the values of -R max which 
corresponds to the radius of neutron star with maximum mass. It is clear from the table 
that for smaller ( the value of Ria varies with Ar and is independent of X^y- This is due 
to the fact that for smaller (, central density for neutron star with mass 1.4M is lower or 
almost equal to the threshold density for hyperons (see also Fig. [7]). But, as ( increases, the 
central density becomes larger than the threshold densities for various hyperons, thus, R±a 
depends on Ar as well as X^y- In Fig. [10] we plot the variations of i? max and Z max versus 
the maximum neutron star mass. We see that the correlations in the values of M max and 
-Rmax are stronger than the ones observed in the case of M max and Ria- The spread in the 
values of -R max is only 0.2 ± 0.1 km for a fixed value of M max . The values of -R max do not 
depend strongly on the choice of Ar as can be seen from Table VI. The horizontal line in the 
lower panel corresponds to the measured value of the redshift, Z = 0.35, for the neutron star 



EXO 0748-676 



53| | . For Z = 0.35, we find that the M max is ~ 1.8M and the corresponding 



radius is ~ 12 km. These values for neutron star masses and the corresponding radii are in 
reasonable agreement with the best suggested value of the mass 1.8M Q and radius 11.5 km 
corresponding to Z=0.35 \&\ - 

To this end, we would like to mention that the calculations are repeated for an attractive 
£ — N potentials by assuming U^ N) = -30 MeV in Eq. (Ep. We find that with this choice 
our results for the variations in the M max and Ria do not get affected. However, 
the threshold density for E~ hyperon becomes lowest and S~ hyperon does not appear even 
for the maximum neutron star mass. It must be pointed out that the tensor coupling of 
cu-meson to the hyperons, not considered in the present work, could increase the value of 
M max by about 0.1M Q [55j. We also remark that the effects due to the exchange and the 
correlations are not considered explicitly. But, they are taken into account at least partly 
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through the non-linear self and mixed interactions of the mesons [29|, |30j • The Eqs. (fTTj) - 
ffT6]) can be interpreted as the Khon-Sham equations in relativistic case and in this sense they 
include effects beyond the Hartree approach through the non-linear couplings. However, a 



more accurate treatment of the exchange and correlation effects should be pursued 
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VI. ROTATING NEUTRON STARS 

The properties of neutron stars can get significantly affected in the presence of rotation. 
The effects of rotation on the neutron star properties are pronounced when the frequency of 
rotation is close to its Keplerian limit. Earlier studies indicate that the Keplerian fre quen cy 
is ~ 1000 Hz for the neutron stars with mass around 1M 



581 ] . Only very recently lOj, a 
neutron star rotating at 1122 Hz is discovered in the X-ray transient XTE J1739-285. In 
this section we shall discuss our results for the rotating neutron stars obtained using the 
extended FTRMF model. These results are obtained by solving the Einstein equations for 
stationary axi-symmetric spacetime. The numerical computations are performed using the 



code written by Stergioulas [59j. 



In Fig. [TT] we plot the neutron star mass versus the circumferential equatorial radius R eq 
for the Keplerian sequences obtained using EOSs for the LY, UY, L0 and U0 parameteri- 
zations of our model. The maximum mass of the neutron stars vary between 1.7 — 2.5M 
and 2.0 — 3.0M Q for the cases with and without the hyperons respectively. The values of 
i?g q 4 lie in the range of 18.4 — 20.0 km irrespective of whether or not hyperonic degrees of 
freedom are included. Because, the central density for the canonical mass in the presence 
of rotation becomes lower than the threshold densities for the hyperons. The Keplerian fre- 
quencies at maximum neutron star mass for various cases shown in Fig. [11] lie in the range 
of 1320 — 1560 Hz . This means that all the EOSs obtained in the present work can yield 
neutron stars rotating at 1122 Hz. In Fig. [12] we plot the mass and the corresponding values 
for i? e q for the neutron star rotating at 1122 Hz. The lower and upper bounds on the radii 
R eq are determined by the setting-in of the axi-symmetric perturbation and mass-shedding 
instabilities, respectively [60j. The maximum values of R eq are well fitted by |60j 

iC q ax = 13.87 (-] km (35) 
which can be obtained using z/& = 1122 Hz and i] = 1 in Eq. (I32I) . In Table IVTlj we give 
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the minium and maximum values for the R eq and the corresponding neutron star mass for 
the various cases plotted in Fig. [121 We get R™ n = 12.1 - 13.8 km and R™ x = 16.5 - 18.7 
km. The values of M(R™ n ) and M(R™ X ) are in the range of 1.6 - 2.7M and 1.7 - 2.6M , 
respectively. The absolute difference between the M(i?™ m ) and M(i?™ ax ) which gives the 
variations in the neutron star mass for a given EOS is at most O.2M . We also find that the 
baryonic mass for the neutron stars rotating with 1122 Hz for all the cases considered here 
are larger than the maximum baryonic mass for the corresponding non-rotating sequences. 
This suggests that the recently discovered fastest rotating neutron star rotating with 1122 
Hz is supra massive. We also list in Table IVHI the values for the r po i e /r eq known as the 
flattening parameter and T j \ W | where T the kinetic energy and W the gravitational 
energy. 

VII. CONCLUSIONS 

We have used the extended FTRMF model to obtain a new set of EOSs with the different 
high density behaviour. These EOSs are then employed to study the variations in the 
properties of non-rotating and rotating neutron stars. The high density behaviour of the EOS 
which is not yet well constrained is varied by choosing the different values of the cu-meson self- 
coupling and the couplings of tu-meson to the various hyperons in our model. The different 
values for these couplings are so chosen that they span the entire range as often considered in 
the earlier works. The remaining parameters of the models are calibrated to yield reasonable 
fit to the bulk nuclear observables, nuclear matter incompressibility coefficient and hyperon- 
nucleon potential depths. The properties of finite nuclei and nuclear matter associated 
with each of the parameterizations used for obtaining EOSs can be summarized as follows. 
The rms errors for the total binding energies and charge radii calculated for the nuclei 
considered in the fits are 1.5 — 1.8 MeV and 0.025 — 0.040 fm. The binding energy per 
nucleon is 16.11 ± 0.04 MeV, saturation density is 0.148 ± 0.003 fm -3 and nuclear matter 
incompressibility coefficient is 230.24 ± 9.80 MeV. 

The values of M max for the non-rotating neutron stars composed of nucleons and hyperons 
in (3 equilibrium can vary between 1.4 — 2.1M . The radius -R1.4 for neutron star can vary 
in the range of 11.3 — 14.1 km. The values of -R1.4 narrow down to only 12.8 — 14.1 km if one 
considers the EOSs for which M max is larger than 1.6M as it is the highest mass measured 
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for PSR J0751+1807 with 95% confidence limit. This result is in reasonable agreement with 
-R1.4 = 14. 81^6 km as deduced very recently by adequately fitting the high quality X-ray 
spectrum from the neutron star X7 in the globular cluster 47 Tucanae [52]. We also note 
strong correlations between the values of M max and -R1.4. The values of the redshift for the 
neutron stars with canonical and maximum masses are also calculated. The redshift for the 
neutron star with canonical mass obtained for different EOSs varies in between 0.19 — 0.25. 
The maximum value for the redshift is 0.41 which corresponds to the maximum neutron star 
mass of 2.1M Q . For the measured value of redshift equal to 0.35, we find that the neutron 
star mass is ~ 1.8M and the corresponding radius is ~ 12 km. These values for neutron 
star masses and the corresponding radii are in reasonable agreement with the best suggested 
value of the mass 1.8M Q and radius 11.5 km for Z=0.35 [54(. For the sake of comparison we 
have presented our results obtained without the inclusions of the hyperons. In this case the 
M max and -R1.4 lie in the range of 1.7 — 2AM & and 12.5 — 14.1 km, respectively. 

We use our EOSs to compute the properties of the rotating neutron stars. In particular, 
we studied the mass and the circumferential equatorial radius for the neutron star rotating 
at 1122 Hz as recently observed [lOj. Our results for different EOSs indicate that the 
mass for such a star can lie within 1.6 — 2.7M . The minimum values for the circumferential 
equatorial radius determined by the onset of the instability with respect to the axi-symmetric 
perturbation are found to vary in the range of 12.1 — 13.8 km. The maximum values for the 
circumferential equatorial radius obtained by the mass-shedding limit vary within 16.5 — 18.7 
km. Looking into the results for the baryonic mass we find that the neutron star rotating 
at 1122 Hz are supra massive for our EOSs. 
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FIG. 1: (Color online) Variations of X a y with X^y for A, E and H hyperons. The values of X a y 
for a given value of X^y are calculated using Eqs. ([241 and |26|) . 



FIG. 2: (Color online) Variations of the symmetry energy coefficient J (upper panel) and its 
linear density dependence L (lower panel) as a function of Ar for different parameterizations with 
C = 0.0,0.03 and 0.06 . 

FIG. 3: (Color online) Variations of the rms errors in the total binding energies (upper panel) and 
charge rms radii (lower panel) as a function of Ar for different parameterizations with £ = 0.0, 0.03 
and 0.06. 

FIG. 4: (Color online) The EOSs for pure neutron matter (upper panel) and symmetric nuclear 
matter (lower panel). The solid and dashed curves correspond to Ar = 0.16 and 0.28 fm, respec- 



tively. The shaded regions represent the experimental data taken from Ref. 



6lfi. 



FIG. 5: (Color online) The chemical potentials for neutron (upper panel) and electron (lower panel) 
as a function of density. 

FIG. 6: (Color online) Variation of neutron star mass as a function of its radius R for selected EOSs. 
These EOSs are obtained using the parameter sets LY (UY) correspond to combinations of Ar = 
0.16(0.28), C = 0.06(0.0) and X^y = 0.5(0.8). The parameter sets L0 and U0 are analogous to 
LY and UY respectively, but, with no hyperons. The various constraints as indicated by causality, 
rotation, and AI / 1 = 0.014 are discussed in the text. 

FIG. 7: (Color online) The threshold density for various hyperons and central densities for neutron 
star with the canonical mass and maximum mass obtained for the EOSs corresponding to the 
parameter sets LY and UY. 

FIG. 8: (Color online) Particle fractions as a function of radial coordinate of the neutron star 
obtained at maximum mass for LY (upper panel) and UY (lower panel) parameterizations. The 
curves labeled as "n/3" should be multiplied by three to get the actual neutron fractions. 



22 



FIG. 9: (Color online) Variations of radius (-R1.4) and redshift (Z1.4) for neutron star with the 
canonical mass as a function of maximum neutron star mass obtained for the EOSs corresponding 
to the all different parameterizations of the extended FTRMF model as considered. The vertical 
line at M max = 1.6M & in the upper panel corresponds to the mass of the PSR J0751+1807 measured 
with 95% confidence limit. 



FIG. 10: (Color online) Variations of radius (R max ) arid redshift (Z max ) for neutron star with the 
maximum mass as a function of maximum neutron star mass obtained for the EOSs corresponding 
to the all different parameterizations of the extended FTRMF model as considered. The horizontal 
line in the lower panel corresponds to the measured value of the redshift, Z = 0.35, for the neutron 
star EXO 0748-676 [53]. 



FIG. 11: (Color online) The relationship between mass M and the circumferential equatorial radius 
R eq for Keplerian sequences for different EOSs obtained within the extended FTRMF model. 

FIG. 12: (Color online) The mass M verses the circumferential equatorial radius R eq for the 
neutron stars rotating at 1122 Hz for selected EOSs obtained within the extended FTRMF model. 
The minimum values for the radius indicated by open circles are determined by the setting-in of 
the instability with respect to axi-symmetric perturbations. The maximum values for the radius 
indicated by open squares are determined by the mass-shedding instability. The values of maximum 
radius are well fitted by solid curve obtained using Eq. ([35|h 
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TABLE I: New coupling strength parameters for the Lagrangian of the extended FTRMF model as given in Eq.([T]). The seven different 
parameter sets correspond to the different values of the neutron skin-thickness Ar for the 208 Pb nucleus used in the fit. The value of w-meson 
self-coupling ( is equal to 0.0 for all these parameterizations. The values of Ar are in fm, the parameters k, oFf, and aFJ are in fm _1 and mo- 
are in MeV. The masses for other mesons are taken to be m w = 782.5 MeV, m p = 763 MeV, m* = 975 MeV and = 1020 MeV. For the 
masses of nucleons and hyperons we use M N = 939 MeV, M A = 1116 MeV, M E = 1193 MeV and M H = 1313 MeV. The values of k, A, aj, 
ax', a~2, a' 2 , and a' 3 are multiplied with 10 2 . 
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TABLE II: Same as Table HI but, with w-meson self coupling £ = 0.03. 
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0.00052 


0.04411 


0.01669 


a.2 


2.89393 


3.06821 


3.18222 


2.77747 


1.18745 


0.68168 


0.47146 


a' 2 


1.59659 


1.16255 


0.47540 


0.22126 


1.27574 


0.54787 


0.52816 


a' 3 


1.52088 


1.35981 


0.97721 


0.45581 


0.28975 


0.35906 


0.32358 


m a 


506.50582 


500.51106 


499.52635 


497.20745 


499.12460 


495.18211 


494.93882 



TABLE III: Same as Table HI but, with w-meson self coupling £ = 0.06. 



Ar 


0.16 


0.18 


0.20 


0.22 


0.24 


0.26 


0.28 


9aN 


11.05170 


11.02412 


10.95765 


11.01908 


10.91944 


11.10806 


11.03151 


9u)N 


14.65579 


14.66595 


14.59582 


14.77458 


14.64700 


15.19792 


15.01572 


9pN 


14.98725 


14.52186 


13.41111 


11.94837 


10.71055 


10.08835 


10.00666 


~K 


0.66576 


0.69497 


0.76852 


0.78002 


0.90221 


1.13349 


0.80797 


A 


2.46427 


2.44874 


2.41259 


2.47238 


2.33265 


2.51229 


2.41320 


Oil 


0.00601 


0.00449 


0.00409 


0.01469 


0.03499 


0.14153 


0.02073 


oil 


0.00203 


0.00526 


0.01079 


0.01559 


0.00230 


0.00085 


0.01109 


a.2 


2.86236 


2.58355 


2.66308 


2.02292 


1.24695 


1.18538 


0.55325 


a' 2 


1.55176 


1.56881 


1.30876 


0.90169 


0.77919 


0.27422 


0.16326 


a' 3 


1.55307 


1.58487 


0.84916 


0.96305 


0.74863 


0.40699 


0.72768 


m a 


503.43838 


501.37038 


499.38134 


497.27203 


495.82388 


490.83495 


490.68907 



TABLE IV: The values of central baryon density p c , mass M, radius R , radiation radius Roo, 
binding energy Lbind an d redshift Z for non-rotating neutron stars with maximum mass calculated 
for the EOSs obtained using LY, UY, LO and UO parameterizations. The parameter sets LY (LO) 
and UY (UO) yield softest and stiffest EOS with (without) hyperons in comparison to all other 
parameterizations obtained in Sec. IIIII 





LY 


UY 


L0 


U0 


Pc (fm 3 ) 


1.05 


0.84 


1.12 


0.79 


M(M ) 


1.4 


2.1 


1.7 


2.4 


R(km) 


11.3 


12.0 


10.9 


12.2 


.Roo (km) 


14.2 


17.3 


14.9 


18.9 


Ebind (10 53 ergs) 


1.36 


3.80 


2.76 


6.49 


Z 


0.25 


0.41 


0.37 


0.57 



TABLE V: Same as Table IV, but, for the non-rotating neutron stars with canonical mass. 





LY 


UY 


L0 


U0 


p c (fm 3 ) 


1.05 


0.32 


0.50 


0.32 


fl(km) 


11.3 


14.1 


12.5 


14.1 


i?oo(km) 


14.2 


16.8 


15.3 


16.8 


L bi „d(10 53 ergs) 


1.36 


1.10 


1.37 


1.10 


Z 


0.25 


0.19 


0.22 


0.19 
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TABLE VI: Values of the M max , R 1A and the radius i? max for the neutron star with maximum 
mass obtained for the EOSs corresponding to the selected combinations of Ar, £ and X u y 









C = o.o 






C = 0.06 




X w y 


Ar 


M max 






M max 


RlA 






(fm) 


(M Q ) 


(km) 


(km) 


(M ) 


(km) 


(km) 


0.50 


0.16 


1.8 


13.4 


12.0 


1.4 


11.3 


11.3 




0.28 


1.8 


14.1 


12.2 


1.4 


11.6 


11.6 


0.80 


0.16 


2.1 


13.4 


12.0 


1.5 


12.3 


11.0 




0.28 


2.1 


14.1 


12.1 


1.5 


13.0 


11.3 



TABLE VII: The properties of neutron star rotating with 1122 Hz for different EOSs calculated 
within FTRMF model. 



EOS 


M{R™ n ) 






T/\W\ 




(MO) 


(km) 






L0 


1.908 


12.12 


0.804 


0.054 


UO 


2.721 


13.46 


0.815 


0.063 


LY 


1.624 


13.85 


0.692 


0.076 


UY 


2.266 


12.99 


0.799 


0.059 


EOS 


M{Rf™) 


pmax 


r pole/ r eq 


T/\W\ 




(MO) 


(km) 






L0 


1.909 


17.14 


0.566 


0.109 


UO 


2.556 


18.73 


0.556 


0.127 


LY 


1.694 


16.49 


0.575 


0.097 


UY 


2.360 


18.35 


0.559 


0.118 
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